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: ABSTRACT 

It has long been speculated that the observed periodic radial velocity pattern for the K giant 
Pollux might be explained in terms of an orbiting planetary companion. We have collected 80 
high-resolution spectra for Pollux at Lick Observatory yielding precise radial velocities with a 
mean error of 3.8 m/s, providing the most comprehensive and precise data set available for this 
star. Our data confirm the periodicity previously seen in the radial velocities. We derive a period 
of 589.7±3.5 days and, assuming a primary mass of 1.86 Mq, a minimum companion mass of 
2.9±0.3Mjup, consistent with earlier determinations. No evidence for any periodicities is visible 
in our analysis of the shapes of the spectral lines via the bisector method, so that we conclude 
that evidence is accumulating and compelling for a planet around Pollux. However, some last 
doubt remains about this interpretation, because non-radial pulsations which might be present in 
giant stars could in principle also explain the observed radial velocities, while the accompanying 
bisector variations might be too small to be detectable with current data. 

Subject headings: stars: individual (Pollux) — planetary systems — techniques: radial velocities — line: 
profiles — stars: oscillations 
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1. Introduction 

Pollux (/3 Gem, HR 2990, HD 62509, HIP 37826) 
is one of the brightest stars in the sky (V=1.16 mag) 
and has been observed extensively in the past. 
Fundamental parameters from a detailed model 
atmosphere analysis of the spectrum have e.g. 
been provided by Drake & Smith (1991), and it is 
usually classified as KOIIIb star (Keenan & Mc- 
Neil 1989). The parallax determined by Hipparcos 
results in a distance of 10.3±0.1 pc. In the Hip- 
parcos Catalogue Pollux was flagged as a possible 
micro-variable with a photometric amplitude of 
less than 0.03 mag (but no obvious periodicity), 
as well as a possibly non-single star, maybe be- 
cause of slightly difi^erent astrometric solutions 
from the two different data reduction consortia. 

Walker et al. (1989) were the first to report sig- 
nificant radial velocity (RV) variations for Pollux, 
with a standard deviation of 26 m/s around the 
mean from RV measurements spread over about 
five years. Though they noted that based on a 
periodogram analysis significant periodicity was 
present in the data, they did not quote any pe- 
riod. Only after having monitored Pollux exten- 
sively over twelve years with a typical RV precision 
of 10-20 m/s, Larson et al. (1993) published a RV 
period of 584.65±3.3 days and discussed possible 
reasons for the observed periodicity. Possible ex- 
planations include an orbiting planetary compan- 
ion or rotational modulation of surface features. 
The latter hypothesis was supported by a slight in- 
dication in the data for a periodicity in the equiv- 
alent width index data of the 8662A (Ca il) line 
with about the same period as found in the radial 
velocities, but with a very low amplitude and only 
a marginal statistical significance. 

Finally, Hatzes & Cochran (1993) presented 
again strong evidence for a periodicity in the radial 
velocities with a period of 558 days. The spectra 
were taken over a period of 3.5 years, and the typ- 
ical RV accuracy was 20 m/s. The RV variations 
were consistent in amplitude and phase with the 
older data by Walker et al. (1989). 

Here we present again precise radial velocity 
measurements of Pollux which leave no doubt 
about a periodicity, determined from our data to 
589.7 days. This RV set is the most comprehensive 
and precise one taken so far for this star, spanning 
almost six years. From the first measurements of 



Walker in 1981 to the latest ones by us in 2006, 
this adds up to 25 years of RV monitoring for Pol- 
lux, with no evidence for a change in phase or am- 
plitude of the almost sinusoidal variations. Along 
with no detectable variations in the spectral line 
shapes, our data set thus lends further evidence 
for the companion hypothesis. 

In Section 2, we describe our observations which 
are part of a larger program of monitoring giant 
stars for periodic RV changes and present our or- 
bital fit to the RV data. In Section 3, we ana- 
lyze the spectral line shapes with the help of the 
bisector method, and in Section 4 we present a 
discussion and our final conclusions. 

2. Observations 

The observations were carried out as part of a 
larger program measuring precise radial velocities 
of several hundred G and K giant stars at Lick 
Observatory. The early objectives of this program 
have been described in Frink et al. (2001), and the 
first substellar companion from the survey (around 
the K giant l Dra) was announced in Frink et al. 
(2002). Hekker et al. (2006) characterize the giant 
stars with the most stable radial velocities in our 
survey, and Mitchell et al. (2007) present evidence 
for four more K giant stars harboring one or more 
substellar companions. 

As part of this ongoing program, we have 
obtained 80 spectra for Pollux covering about 
5.5 years. All observations were taken using the 
Hamilton High Resolution Echelle Spectrograph 
(R«60 000 at 6000 A) at Lick Observatory, at- 
tached to the 0.6 m Coude Auxiliary Telescope 
(CAT). Typical exposure times were 1.5 minutes 
for Pollux, yielding a S/N of up to 200. However, 
some observations were taken with cloud cover, 
where exposure times can be considerably longer 
and the S/N somewhat smaller (around 150). The 
individual radial velocities, obtained as described 
in Butler et al. (1996), are listed in Table 2, along 
with their formal errors. Figure 1 shows the mea- 
surements together with a Keplerian fit to the 
data; the orbital elements are listed in Table 1. 
There is no doubt about the clear periodicity in 
the data, as has already been observed by Walker 
et al. (1989), Larson et al. (1993) and Hatzes & 
Cochran (1993). 

Nevertheless, there is some additional scatter at 
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Table 1 

Fitted and derived orbital parameters 





V aJuc 


Period [days] 


589.7 ±3.5 


To [JD-2 450000] 


2337.9t^^ 


Eccentricity 


0.06 ± 0.04 


w [deg] 


277 ±8 


/(m) [IO-^Mq] 


6.2 ±0.6 


7712 smi [Mjup]'' 


2.9 ±0.1 


Semi-major axis [AU] 


1.69 ±0.03 


RV semi-amplitude [ms~^] 


46.9 ± 1.5 


Reduced 


6.3 


rms scatter around fit [ms~^] 


9.0 



'^Thc companion mass error docs not include 
the uncertainty in the stellar mass. 



a level of 9 m/s present in the data, which is larger 
than expected based on the formal measurement 
errors. It is possible that this stems from radial 
pulsations, but with the theoretical period for the 
fundamental mode shorter than 1 day (Hatzes & 
Cochran 1993) our sampling is inadequate to pro- 
vide any further constraints. Solar-like oscillations 
in late G giants have been found by Frandsen et 
al. (2002) and de Kidder et al. (2006), and it is 
not unreasonable to assume that similar oscilla- 
tions in Pollux are responsible for the excess jit- 
ter. The amplitude however is not large enough to 
affect the derivation of the orbital parameters of 
the putative companion. Furthermore, there are 
no indications for any additional periodicities, as 
can be seen in the lower part of Fig. 1 where the 
orbital fit has been subtracted from the data, and 
the corresponding Lomb-Scargle periodogram in 
the lower part of Fig. 2. 

Our period for the orbital fit is 589.7±3.5 days 
and compares very well with the value of Larson 
et al. (1993) (584.65±3.3 days), while it is some- 
what larger than the period derived by Hatzes & 
Cochran (1993) (554±8 days). All previous mea- 
surements arc consistent in amplitude and phase 
with our result; Fig. 3 shows a phased plot of our 
measurements along with the earlier ones by Lar- 
son et al. (1993) and Hatzes & Cochran (1993). 
(Note that we measure only relative radial veloc- 



ities with an arbitrary zero point, so that a ver- 
tical shift had to be applied before plotting them 
along with the other data sets. This vertical shift 
is a free parameter in our fit; its formal error is 
1.1 m/s and should be negligible for the compari- 
son.) A combined fit to all RV data produces very 
similar orbital parameters as a fit to our RV data 
alone; the period from the combined fit is larger 
by 1.6 days as compared to the period quoted in 
Table 1 based on our RV data alone, very well 
within the formal error. 

The minimum companion mass derived from 
our orbital fit is 2.9 Mjup assiiming a stellar mass 
for Pollux of 1.86 M©. The stellar mass was de- 
rived from the location of the star in the color- 
magnitude diagram as determined from Hippar- 
cos data and compared to the evolutionary tracks 
from Girardi et al. (2000). Solar metallicity was 
assumed for the comparison, which is a good ap- 
proximation (Drake & Smith 1991). AUende Pri- 
eto & Lambert (1999) derive a value of 1.7 Mq for 
the mass of Pollux with a very similar method as 
described above, and Drake & Smith (1991) also 
derive a mass of 1.7 Mq from a detailed model at- 
mosphere analysis. Using this value for the stellar 
mass would yield a minimum mass of 2.7 Mj^p. 
The error on the minimum companion mass due 
to the error in the knowledge of the primary mass 
is thus about 0.2 Mj^p, so that the total error (in- 
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eluding the formal error derived from the orbital 
fit, see Table 1) amounts to about 0.3 Mj^p. 

Together, the radial velocities cover about 
25 years, and the variations have been rather 
consistent over that time. Nevertheless, though 
we consider the interpretation of the observed 
RV changes as the result of an orbiting compan- 
ion likely, it is possible that another mechanism 
might cause the observed RV variations. In Chap- 
ter 3, we take a closer look at the spectral line 
shapes which might provide further hints at the 
underlying mechanism. 




500 1000 1500 2000 2500 

JD-2451545.0 [days] 

Fig. 1. — Top Panel: Radial velocities measured at Lick Observatory, along with error bars, covering about 
5.5 years from September 2000 to February 2006. The best fit Keplerian is overplotted, with a period of 
589.7 days. Bottom Panel: Radial velocity residuals after the best fit Keplerian has been subtracted. The 
remaining radial velocity scatter has a standard deviation of 9 ni/s, and no systematics are visible in the 
residuals, neither by eye nor in the periodogram of the residuals in the lower part of Fig. 2. 
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Fig. 2. — Top Panel: periodogram of the measured radial velocities. The highly significant peak occurs 
at 588 days; a Kepler fit to the data reveals a best fit period of 589.7 days. The next significant peak to 
the right, at about log P[days]=3.07, corresponds to twice the value of the most significant period. The 
numbers at the right indicate the false alarm probabilities of the labeled lines; a highly significant peak in 
the periodogram would sit clearly above the highest line, indicating a false alarm probability of less than 
0.1%. Bottom panel: same as above, but with the Kepler fit corresponding to 589.7 days removed from the 
radial velocities. No significant peak is left in the periodogram of the radial velocity residuals. 
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Fig. 3. — Radial velocities from Larson et al. (1993) (open triangles) and Hatzes & Cochran (1993) (asterisks) 
phased to the period determined from our own measurements (filled circles). Our orbital fit is also shown 
(solid line). No changes in period or amplitude of the periodic RV signal is apparent from that plot. For 
clarity, no individual error bars have been included in the plot. The data in Larson et al. (1993) consist of 
two separate datasets; data from CFHT have mean internal errors of 12 m/s, data from DAO 27 m/s. The 
majority of the datapoints in Hatzes & Cochran (1993) has mean internal errors of around 20 m/s, while 
a few datapoints have mean internal errors of 7 m/s. Note that not all of the measurements in Hatzes & 
Cochran (1993) have been used; the second set of 13 measurements in their Table IC is identical to the 
first 13 measurements and thus obviously wrong, so it has been omitted. Errors for the DAO dataset from 
Larson et al. (1993) are rather inhomogeneous and reach up to around 80 m/s for some measurements, which 
explains the few data points which seem to be outliers in this graph. 
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3. Line Shape Analysis 

In order to investigate whether the observed RV 
variations are caused by a shift of the spectral Unes 
as a whole (as expected in the presence of a com- 
panion) or by a change in the symmetry of the 
spectral lines giving rise to a net change in RV 
(as expected in the presence of pulsations), bisec- 
tors of the cross correlation profile have been an- 
alyzed. We used all spectral lines between about 
6540 A and 9590 A from 29 spectral orders and ob- 
tained an average line profile by correlation with 
a synthetic template which was obtained from 
the VALD database (Kupka et al. 1999, available 
at http:/ /ams. astro. univic.ac.at/vald/), matching 
the effective temperature and surface gravity of 
Pollux. The spectral range from about 5000 A to 
5800 A could not be used because it is affected by 
iodine lines, and the spectral range around 6300 A 
was avoided because it is dominated by strong tel- 
luric oxygen lines (which give rise to a spurious 
one year period when included in the analysis). 
Otherwise, as many lines as possible from a con- 
tinuous spectral range were used, since it is known 
(Gray 1983, 1984) that different lines in the same 
star can display different bisector behavior. 

For the cross correlation with the synthetic 

template, our individual spectra had to be wave- 
length-calibrated, for which we used the Thorium- 
Argon exposure which was taken closest in time to 
each spectrum, either from the beginning or end 
of the night. Altogether, about 1110 spectral lines 
with theoretical depths between 0.1 and 0.9 were 
used for the cross correlation. Of course many 
blended lines are included in the cross correlation, 
but by using many lines we are confident that the 
effects average out. Furthermore, as long as the 
same lines are used for all observations, and only 
variations in the shape are of interest, blends by 
stellar lines do not affect the final result. 

After having obtained the cross correlation pro- 
files, we determined bisectors by stepping down 
the blue side of the profile, linearly interpolating 
the line depth at the red side for the same flux 
level as observed at the blue side (in the center 
of the profile, a parabolic fit was used instead of 
the linear one), and derived the midpoints between 
the velocities on the blue and red side of the pro- 
file. The connection of midpoints determined in 
this way is the bisector; see Povich et al. (2001) 



for more details on the method. 

In order to analyze possible variations in the 
bisector, two quantities are defined which charac- 
terize its shape: the velocity span and the velocity 
displacement. While the bisector span is the dif- 
ference between the width of the bisector at two 
different fiux levels (30% and 75% were used here), 
the velocity displacement is the average width of 
the bisector at three different flux levels (30%, 60% 
and 75% were used). 

Periodograms of both quantities are shown in 
Fig. 4. No significant periods whatsoever are 
present in these periodograms; all trial periods 
have extremely small significance levels. In par- 
ticular, no peak is present at the RV period, so 
that we conclude that there is no evidence for any 
variations in the shapes of the spectral lines in our 
Pollux spectra. 

This finding is consistent with the one by 
Hatzes & Cochran (1998), who also analyzed the 
bisectors of two different spectral lines in their 
Pollux data without discovering any periodicities 
similar to the RV period (they used 554 days). 
However, they caution that low-order non-radial 
pulsation modes, which might be able to account 
for the observed RV variations, might produce 
changes in the bisector velocity span of only 5 to 
20 m/s. While we estimate the errors in a sin- 
gle bisector velocity span to lie between 50 and 
100 m/s in our analysis (50 m/s in the analysis of 
Hatzes & Cochran (1998)), the sensitivity to a real 
periodic signal in the data increases if one uses a 
large number of spectra for the periodogram, to 
about the relevant level. However, we conclude 
that non-radial pulsation cannot be completely 
ruled out with the observational material or anal- 
ysis methods currently available. 
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Fig. 4. — Periodograms of the derived bisectors. Top panel: velocity span, bottom panel: velocity displace- 
ment. There are no significant periodicities visible; all trial periods have an extremely small probability of 
being real. The period present in the RV data is indicated. 



9 



4. Conclusions 

The presently available data on Pollux are all 
compatible with an orbiting substcUar companion 
around this bright and nearby star. The minimum 
mass of the companion is 2.9±0.3 Mjup, which 
makes it most likely a planet and not a brown 
dwarf. From the non-detection of the companion 
in the Hipparcos data, one can derive an upper 
limit on the companion mass, since it would have 
been detected by Hipparcos if it had been massive 
enough. The lower limit on the inclination which 
we derive from Hipparcos is about 5°, which trans- 
lates into an upper mass limit of 33 Mjup, con- 
straining the companion to be of substcUar na- 
ture (if the companion interpretation of the RV 
pattern is correct). Its period is 589.7±3.5 days, 
and it orbits at a distance of 1.69±0.03 AU from 
the star in an almost circular orbit. The RV pat- 
tern has been stable over the last 25 years. An 
analysis of the spectral line shapes shows no ev- 
idence for any changes with the RV period nor 
any other periods, which one might expect in the 
presence of pulsations. Hipparcos has picked up 
excess scatter in the astrometric standard solu- 
tion (without companion), but it is unlikely that 
this is the signature of the planetary companion, 
since its expected minimum astrometric signature 
is only 50 /xas. 

Detailed theoretical predictions of the expected 
amount of bisector asymmetry in the presence of 
non-radial g- or r-mode pulsations in giant stars 
are currently not available. However, the numeri- 
cal simulations conducted by Hatzes (1996) show 
that it is in principle possible to explain the ob- 
served RV variations in Pollux by low-order non- 
radial pulsation modes, while the accompanying 
bisector variations would be too small to be de- 
tected with current techniques (Hatzes & Cochran 
1998). This is the main reason why some last 
doubt remains about the interpretation of the RV 
variations in terms of an orbiting companion, even 
if it would be a bit surprising to find only one 
single long-period pulsation mode. Note that the 
small eccentricity which we find (0.06±0.04) is 
barely significantly different from zero, so that it 
cannot be used to rule out pulsations as the reason 
for the observed RV variations as has been done 
for highly eccentric giant star orbits (Frink et al. 
2002). 



In contrast to that, rotational modulation of 
starspots can be excluded as the reason for the 
observed RV changes, since otherwise some photo- 
metric variability larger than the micro- variability 
actually seen should have been detected by Hip- 
parcos. Also, it would be difficult to explain how a 
single or several starspots could produce RV varia- 
tions which are so close to sinusoidal over the rota- 
tion period as well as stable over the last 25 years. 

We conclude that while evidence is accumulat- 
ing and compelling for an orbiting planet around 
Pollux, the final confirmation has to await a the- 
oretical prediction of the amount of spectral line 
asymmetry in the presence of non-radial g- or r- 
mode pulsations in giant stars, much increased 
sensitivity in bisector analyses or photometry or 
the detection of the companion with independent 
techniques like e.g. precise astrometry. 

After this paper was first submitted we learned 
of the similar paper by Hatzes ct al. (2006). Hatzes 
et al. (2006) present 55 new radial velocity mea- 
surements of Pollux with mean internal errors be- 
tween 11 and 17 m/s taken between 1998 and 2006 
and analyze them together with the older datasets 
by Larson et al. (1993) and Hatzes & Cochran 
(1993) also used in the present paper. Their or- 
bital elements are in excellent agreement with the 
ones derived in the present paper, and both papers 
arrive at the same conclusions regarding the inter- 
pretation of the observed RV periodicity in terms 
of an orbiting planetary companion. 

We kindly thank the staff at Lick Observatory 
for their outstanding dedication and support, as 
well as our referee, Gordon Walker, for helpful 
comments on the manuscript. 

This research has made use of the Vienna 
Atomic Line database (VALD) located at 
http: / / ams.astro.univie.ac.at/vald/. 

REFERENCES 

Allende Prieto, C, & Lambert, D. L. 1999, A&A, 

352, 555 

Butler, R. P., Marcy, G. W., WiUiams, E., Mc- 
Carthy, C, Dosanjh, P., & Vogt, S. S. 1996, 
PASP, 108, 500 

de Ridder, J., Barban, C, Carrier, F., Mazumdar, 



10 



A., Eggcnbcrgcr, P., Acrts, C, Dcruytcr, S., & 
Vanautgacrdcn, J. 2006, A&A, 448, 689 

Drake, J. J., & Smith, G. 1991, MNRAS, 250, 89 

Frandscn, S., et al. 2002, A&A, 394, L5 

Frink, S., Quirrcnbach, A., Fischer, D., Roser, S., 
& Schilbach, E. 2001, PASP, 113, 173 

Frink, S., MitcheU, D. S., Quirrenbach, A., Fis- 
cher, D. A., Marcy, G. W., & Butler, R. P. 2002, 
ApJ, 576, 478 

Girardi, L., Brcssan, A., Bertelli, G., & Chiosi, G. 
2000, A&AS, 141, 371 

Gray, D. F. 1983, PASP, 95, 252 

Gray, D. F. 1984, PASP, 96, 382 

Hatzes, A. P., & Cochran, W. D. 1993, ApJ, 413, 
339 

Hatzes, A. P. 1996, PASP, 108, 839 

Hatzes, A. P., & Cochran, W. D. 1998, ASP 
Gonf. Ser. 134: Brown Dwarfs and Extrasolar 
Planets, 134, 312 

Hatzes, A.P., Cochran, W.D., Endl, M., Guen- 
ther, E.W., Saar, S.H., Walker, G.A.H., Yang, 
S., Hartmann, M., Esposito, M., Paulson, D.B. 
2006, accepted by A&A 

Hckkcr, S., Reffert, S., Quirrcnbach, A., Mitchell, 
D. S., Fischer, D. A., Marcy, G. W., Butler, 
R. P. 2006, accepted by A&A 

Keenan, P. C., & McNeil, R. C. 1989, ApJS, 71, 
245 

Kupka, F., Piskunov, N., Ryabchikova, T. A., 
Stempels, H. C, & Weiss, W. W. 1999, A&AS, 
138, 119 

Larson, A. M., Irwin, A. W., Yang, S. L. S., Good- 
enough, C, Walker, G. A. H., Walker, A. R., & 
Bohlender, D. A. 1993, PASP, 105, 825 

Mitchell, D.S., Reffert, S., Quirrenbach, A., Fis- 
cher, D. A., Marcy, G. W., Butler, R. P. 2007, 
in preparation 

Povich, M. S., Giampapa, M. S., Valenti, J. A., 
Tilleman, T., Barden, S., Deming, D., Liv- 
ingston, W. C., & Pilachowski, C. 2001, AJ, 
121, 1136 



Walker, G. A. H., Yang, S., Campbell, B., & Irwin, 
A. W. 1989, ApJ, 343, L21 



This 2-column preprint was prepared with the AAS lATJ^X 
macros v5.2. 



11 



Table 2 

Measured radial velocities for Pollux 



JD" [days] 


Vrad [ms ^] 






35.4 


3.1 




32, 7 


2.9 


ioo4.U i o 


41.5 


3.1 


1 orto f\£;ri 

loyo.uou 


56.0 


3.8 




57.3 


3.1 


ioyy.csyy 


55.1 


3.2 


1 ciCi-i ono 


50.5 


3.7 


lyzy .oiy 


38.7 


3.5 


1 riQn ooQ 
iyoU-OMO 


42,9 


3.2 


lyoi.Yzo 


34.4 


4.5 


iy4y.77i 


39,4 


5.7 


iyyu.DD ( 


21.2 


3.3 


iyy^.oyo 


32.3 


3.8 


^Ulo.oyy 


13.9 


3.2 


^UoZ.OfSi 


— 6.9 


3.3 


zU4o.OOD 


— 5.9 


3.6 


zU4 / .DO / 


0.3 


3.7 


zU4o.D5y 


3.8 


3.5 


ziYo.uao 


—42.5 


2.9 


zi77.U47 


—48.8 


3.3 


ziyo.u^i 


— 44.7 


3.6 


ziyo.yoo 


— 42.8 


3.2 


zzuo.yoo 


—36.6 


3.3 


z^UD.y7u 


—47.1 


4.1 


oon*? nor 

zzu / .yyo 


—38.4 


3.5 


zzzz. yi / 


—34.3 


3.4 


noHQ nn /I 

zzzo.yz4 


—37.8 


3.1 


Zzoo.ooZ 


—23.9 


3.7 




— 13.9 


4.0 


TOO c; o 1 o 


— 17.2 


3.9 


ZZy t.tilO 


1.5 


4.2 


zou * . *4y 


—0.3 


3.7 


ooen Ton 
zoDZ. t zy 


14,2 


4.3 


2363.740 


14,5 


4.4 


oo o H 'yn 1 
ZOo4.7Ul 


27.9 


4.9 


Zoyo.D f b 


26.2 


3.9 


o CL on n /I Q 
zozy .U4o 


44,8 


3.3 


2542.035 


13,7 


3.8 


2544.004 


14.6 


3.1 


2560.031 


37.6 


5.4 


2561.051 


18.6 


3.2 


2562.023 


29.2 


3.1 


2572.020 


25.0 


3.1 


2590.000 


15.0 


3.4 


2603.974 


8.2 


3.7 


2604.957 


12.3 


3.8 


2605.893 


6.7 


3.5 


2615.882 


2.0 


4.0 


2616.889 


-6.4 


3.7 


2627.889 


-20.9 


3.2 


2665.800 


-25.5 


4.1 


2667.840 


-27.8 


4.7 


2699.669 


-26.2 


4.8 


2765.687 


-35.8 


5.6 


2933.021 


13.1 


3.8 


2935.010 


6.7 


3.0 


2964.043 


42.7 


4.9 


2964.938 


36.3 


4.0 


2966.919 


27.8 


3.4 
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Table 2 — Continued 



JD» [days] 






2985.912 


43.4 


5.8 


3022.829 


40.2 


4.9 


3025.893 


41.6 


4.0 


3089.732 


33.4 


5.3 


3111.687 


37.5 


4.5 


,'S2!)0.027 


-27.1 


3.0 


3323.988 


-46.7 


3.8 


3324.988 


-37.6 


3.9 


3355.851 


-56.1 


3.3 


3357.988 


-55.3 


3.2 


3400.828 


-58.3 


3.3 


3424.734 


-60.8 


3.7 


3442.703 


-49.5 


3.6 


3446.642 


-31.1 


4.0 


3447.665 


-41.5 


4.1 


3492.671 


-15.7 


3.4 


3651.058 


37.7 


3.5 


3703.043 


20.2 


4.3 


3741.009 


22.7 


4.3 


3790.775 


1.9 


3.7 


3792.756 


-0.8 


3.9 



* Julian date - 2 450000 
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